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ABSTRACT 

We investigate the quantitative constraint on the triple-a reaction rate based on stellar evolution 
theory, motivated by the recent significant revision of the rate proposed by nuclear physics calculations. 
Targeted stellar models were computed in order to investigate the impact of that rate in the mass 
range of 0.8 < M/Mq < 25 and in the metallicity range between Z = Q and Z = 0.02. The revised 
rate has a significant impact on the evolution of low- and intermediate-mass stars, while its influence 
on the evolution of massive stars (M > IOMq) is minimal. We find that employing the revised rate 
suppresses helium shell hashes on AGB phase for stars in the initial mass range 0.8 < M/Mq < 6, 
which is contradictory to what is observed. The absence of helium shell flashes is due to the weak 
temperature dependence of the revised triple-a reaction cross section at the temperature involved. In 
our models, it is suggested that the temperature dependence of the cross section should have at least 
i^ > 10 at T = 1 - 1.2 X 10® K where the cross section is proportional to T'''. We also derive the helium 
ignition curve to estimate the maximum cross section to retain the low-mass first red giants. The 
semi-analytically derived ignition curves suggest that the reaction rate should be less than ^ 10~^^ 
cm^ s^^ molc^^ at « 10^'® K, which corresponds to about three orders of magnitude larger than 
that of the NACRE compilation. In an effort to compromise with the revised rates, we calculate and 
analyze models with enhanced CNO cycle reaction rates to increase the maximum luminosity of the 
first giant branch. However, it is impossible to reach the typical RGB tip luminosity even if all the 
reaction rates related to CNO cycles are enhanced by more than ten orders of magnitude. 
Subject headings: stars: evolution — stars: AGB and post AGB — nuclear reactions 



1. INTRODUCTION 

Nuclear reaction rates arc fundamental to the mod- 
eling of stars. Among them, helium-burning reactions 
by the triple-a process play a critical role in evolution 
of red giants. The cross section of the triple-a reaction 
is largely infiuenced by resonances in the energy states 
of '^He, ®Be and ^^C. These resonances were proposed 
bv lOpikI ()195H ) and iSalpeterl (19521); the res onanc e level 
of ^^C was predicted by iHovle et al.1 (|1953[) and iHovld 
(|1954[ ) to be at 7.68 MeV, confirmed ex perimentally with 
a value of 7.68±0.03 MeV (jDunbar et a l. 1953). The con- 
tribution bj;jion;Tesonantreaction was also taken into ac- 
count bv lNomoto et al.l (|1985[ ) to investigate the impact 
on helium ignition in accreting white dwarfs and neu- 
tron stars. Their result significantly increases the triple- 
a reaction rate by more than 20 orders of magnitude at 
r < 3 X 10^ K. However, this drastic change of nuclear 
reaction rates do not affect on the main characteristics of 
stellar evolution because the rate is not so significantly 
modified within the temperature range for helium igni- 
tion in stars {T > 6 x 10^ K). The only one exception 
is the behavio r of the evolutio n of low-mass m etal-free 
stars (,Fu1imo to et al.l [1990; .Suda et al.ll2007aD . With- 



out CNO elements initially, the hydrogen fiash is trig- 
gered by CNO cycles as a consequence of the accumula- 
tion of carbon in the hydrogen-burning domain through 
the triple-a reaction. However, this phenomenon cannot 
currently be observed since the probability of detecting 
this event is significantly small from the viewpoint of 
both timescale and s tatistics. Thus the reaction rate by 
iNomoto et al.l ()1985n has been adopt ed by the later com- 
ilations of nuclear react ion rates (jCaughlan fc Fowler! 



1988: Angulo ct alj Il999f ) in place of previous work, 
which did not include the co ntribution of non- resonant 
reactions (iFowler et al.lll975D. 

RecentlvT IOgata et al.l (|2009[) (OKK, hereafter) pre- 
sented results on the cross section of the triple-a reac- 
tion using a direct calculation of Schrodinger equation. 
They also provide a tremendously larger cross section 
at 10^"® K compared with the previous compilations for 
this nuclear reaction rate. The difference amounts to 
26 orders of magnitude a round 10^ K comp ared with 
the NACRE com pilation (lAngulo et a 1 119991 ) . i.e., the 
reaction rates by INomoto et al.l (11985), an d also with 



the latest reaction rate (jFvnbo et alJ 



2005f ). This re- 
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suit gives rise to a severe inconsistency with the cur- 
rent understandings of the observations. In particular, 
low-mass red giants are no longer produced using the 
revised triple-a reaction rate due to the premature igni- 
tion of helium at the stellar center upon arrival at the 
base of the red giant br anch ([Potter fc Pa xton 20091 . In 
addition to th i s repo rt, iSaruwatari fc Hashimoto (2010 ) , 
iMatsuo et all (|201in . and lPeng fc OttI (I2O1O0 hwestigate 
the impact on accreting C-0 white dwarfs and neutron 
stars. For the accretion onto white dwarfs, it is suggested 
that the energy released during off-center helium detona- 
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tion attained when employing the OKK rate dominates 
over that of the central carbon ignition that is responsible 
for Type la supernovae. On the other hand, it is reported 
that accretion onto neutron stars with the OKK triple-a 
reaction rate results in burst energies from ultracompact 
X-ray b inaries that ar e too l ow to account for those ob- 
served. iMatsuo et al.l ()2011[ ) insist that the OKK rate 
should be reduced by a factor of 10^~^ for lO^K to be con- 
sistent with the observations of X-ray burs ts. Applica- 
tion of the OKK rate to Cepheid models by iMorel et al.l 
(|201Qil also fails to explain the accumulations of obser- 
vations reporting the period change, depletion of lithium 
and reduction of the C/N ratio, although they insist that 
the OKK rate may reduce the mass discrepancy between 
stellar evolution models and pulsation theory. 

Stimulated by this revision of the triple-a reaction 
rate, new calculations of triplc-a reaction rates such 
as a method solving the Faddeev three-body equations 
have been ongoing (S. Ishikawa, private communication). 
It is also proposed that broad resonances other than 
the Hoyle state exist, which may also affect the triple- 
a reaction at low temperat ures (see, e.g . . iKato et afl 
120100 . On the other hand, lOgata et all ([20091) point 
out that the treatment o f non-resonant r eactio ns at low- 
temperatures adopted bv lNomoto et al.l (|1985[ ) is not ap- 
propriate. Therefore, it is desired to explore the reaction 
rates both theoretically and experimentally, and hence, 
it is likely that the cross section for triple-a reaction will 
be modified by future works. As a calibration test for 
the modification of the reaction rate, it is useful to rely 
on stellar evolution models by analyzing the dependence 
of the theoretical models on the reaction rate, especially 
below 3 X 10^ K where the influence is important for the 
evolution of low- and intermediate-mass stars. 

In this paper, we provide quantitative constraints on 
the triple-a reaction cross section based on stellar evolu- 
tion theory. Firstly, we investigate the impact of chang- 
ing the triple-a reaction rate on the evolution of low- 
and intermediate-mass stars during the horizontal branch 
and AGB phase, extending our study to the evolution- 
ary characteristics of massive stars. The possible ef- 
fects of changing the triple-a reaction rate are analyzed 
by comparing our models with the observed extremely 
metal-poor stars in the Galaxy using the SAGA database 
(jSuda et al.l l2008| ). Secondly, we derive the possible 
range of the triple-a reaction rate by estimating the igni- 
tion point of helium burning on the first red giant branch 
(FGB) in a semi-analytic method. We also explore the 
possibility of varying the input physics with the OKK 
triple-a reaction rate adopted. In particular, we discuss 
the change of the reaction rates of CNO cycles. 

The paper is organized as follows. In the next section, 
the stellar evolution models used in this work are de- 
scribed. The impact of changing the triple-a reaction 
rates is also discussed. § E] is devoted to the deriva- 
tion of helium ignition curve and the comparison with 
observations. In § |4l we look for a compromise with 
the OKK triple-a reaction rate by varying other input 
physics. Summary and conclusions follow in § [5] 

2. STELLAR EVOLUTION WITH THE OKK RATE 

The stellar evolut ion code and adopted in put physics 
are the same as in iSuda fc Fuiimotol (I2010D except for 
conductive opacities for which we adopt the calculation 



bv lltoh et al.l (|2008f ). For the nucle ar reaction rates, we 
use both the NACRE compilation (jAngulo et al.lll999D 
and the OKK rate. The mixing length parameter a is set 
at 1.5. We followed the evolution from zero-age main se- 
quence to the early phase of the AGB. Initial c omposition 
is sca led solar abundance based on lAnders fc Grevessd 
(|1989() without any a-element enhancement. We also 
discuss the effect of using the OKK rate on the evolu- 
tion of massive stars and present the results with the 
mode l of 12 and 25M(7i w ith Z = 0.01 using the Geneva 
code (jHirschi et al.l 12004 ). For these models, a = 1.6 is 
adopted. For massive star models, calculations are ter- 
minated at the beginning of carbon burning phase. Here- 
after we refer to the models computed with NACRE rates 
as fiducial models. 

Figure [T] compares evolutionary tracks for various sets 
of mass and metallicities using the NACRE and the OKK 
rate. For models of M > IOMq, the impact of adopting 
the OKK rate is not significant since more massive stars 
ignite helium at higher temperatures and lower densities. 
At the start of core He burning, the nuclear burning rate 
is larger in the OKK rate than in the NACRE rate for a 
12Mq model by a factor of ~ 10"^, while it is only larger 
by a factor of 5 to 10 for a 25Mq model. 

The helium ignition temperature is 1.14x lO'^ K (1.55x 
10^ K with the NACRE rate) and 1.52 x 10^ (1.75 x 
10^) K for 12Mq and 25Mq respectively, and is weakly 
dependent on the triple-a reaction rates adopted for this 
mass range. The OKK rate gives larger luminosity for 
the 12M0 model, placing it higher in the H-R diagram 
because helium burning starts in contraction phase, while 
it occurs on red giant branch for the fiducial model. The 
OKK rate also gives larger final ^^C mass fraction at the 
end of core helium burning. A higher ^^C mass fraction 
can change the size of the convective cores during the 
advanced burning stages and also the final iron core mass. 
For 25M0, such differences are negligible. 

For low- and intermediate-mass stars, the effect of us- 
ing the OKK rate is significa nt even after core helium 
burning. As pointed out by iDotter fc PaxtonI p009) . 
the helium ignition occurs in the center when the star 
is at the base of the FGB on the H-R diagram for 
an initial mass of O.SAfQ, while a helium flash occurs 
off-center for the initial mass of M/Mq < 2 for the 
fiducial models, with some dependence on metall icity 
(see. e.g. lSuda et al.l[2007bl : iSuda fc Fuiimotol [20ll . In 
these low-mass stars, the luminosity at the horizontal 
branch phase is smaller by a factor of 3-60 compared 
with fiducial models because of small helium core masses 
(~ 0.2-0.25A'/q) at the point of helium ignition. Be- 
cause of the low luminosity, the duration of the core 
helium burning can be longer than the fiducial models. 
Since the typical lifetime on horizontal branch is 50- lOO 
Myrs, low-luminosity horizontal branch stars should be 
detected in the observations of globular clusters if the 
OKK rate is adopted. Once the core has exhausted he- 
lium in the center, the star ascends the red giant branch 
for the first time. This can be interpreted as red giants, 
but there is a significant gap between the helium igni- 
tion point and the ascension of the giant branch, which 
has never been detected in the observed color-magnitude 
diagrams. 

Another remarkable difference is that the models with 
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M > 0.8 A/© do not experienee thermal pulsations of he- 
lium shell burning during the AGB phase. Figure [2] shows 
the variations in helium burning luminosity as a funetion 
of helium core mass during AGB evolution. Models are 
taken with the NACRE or OKK rate for the mass range 
of O.8-5M0. In the OKK models, helium shell-burning 
occurs in a radiative zone with a front that advances 
steadily without thermal pulses. This leads to a smooth, 
gradual increase of the carbon-oxygen core mass together 
with a gradual increase of the helium burning luminosity. 
The absence of thermal pulses or shell-flashes is a conse- 
quence of the weaker temperature dependence of triple-a 
cross sections for the OKK rate. Figure [3] compares the 
temperature dependence, v, of triple-a reaction rates, 
defined by i^ = d\og{av)aaa/d\ogT, as well as the cross 
section itself. The value of ly for the OKK rate is smaller 
than ?» 8 for the temperature range for helium shell burn- 
ing (1-1.2 X 10* K) in the present computations. The 
absence of thermal pulses is a serious challenge to obser- 
vations. Without thermal pulses, AGB models do not 
experience any third dredge-up events that are respon- 
sible for the enhancement of carbon and s-process ele- 
ments. Therefore, the mass transfer scenario, which is 
likely to account for the carbon-enhancement observed 
in CH and CE MP stars (includin g CEMP-s stars), is not 
viable (see, e.g. lSuda et al.|[20lll ). 

The operation of the third dredge-up is also sup ported 
by the detection of ^^^Te in AGB stars (iMerrilll [TOSl 
iSmith fc WallersteinI HOSl . Since the half-life of ^^Tc 
is only 2.13 x 10^ years, this element is thought to be 
produced by the s-process and dredge d-up to the sur - 
face during each thermal pulse cycle (iCameronI I1955D . 
In addition to avoiding dredge-up events, we also ex- 
clude the possibility for efficient mass loss, such as su- 
per wind, triggered during the late stage of AGB evolu- 
tion because of the lack of carbon that c an be condensed 
into dust to drive the wind mass loss (llben fc Renzinil 
[1981 [Boron fc WillsonlllMI: iLaeadec fc Ziilstrall2008D . 
Without the efficient mass loss, AGB stars can not 
evolve to white dwarfs, and hence explode as supernovae 
Type 11/2 when the core mass approache s the Chan- 
drasekhar mass limit (llben fc Renzini|[l983D . 

With regards to the subsequent stage of evolution, the 
models with the OKK rate ignite carbon for smaller 
initial mass (M > 6Mq) compared with the mod- 
els with the NACRE rate (M > 8Mq) (see Fig.l of 
ISuda fc Fuiimotoll2010D . This is because the OKK mod- 
els have a larger core mass after the second dredge-up 
(see Fig. [5]) due to the earlier helium ignition during 
the contraction phase after the helium core exceeds the 
Schonberg-Chandrasekhar mass limit. In addition, the 
carbon-oxygen core also becomes massive enough to ig- 
nite carbon in the core before the second dredge-up re- 
duces the helium-rich layer. In our models of 7 and 
9Mq with Z = 0.01, the mass of the carbon-oxygen 
core becomes as large as 1.07 and 1.27Afo when the 
OKK rate is adopted, respectively. These are in con- 
trast with the models with the NACRE rate for which 
we find 0.89 and I.OIM© for 7 and 9Mq models re- 
spectively. Since the carbon-oxygen core is weakly de- 
generate for 7 < M/Mq < 9 model with the OKK 
rate, these models will becom e so-called super- AGB stars 
(jGarcia-Berro fc Ibenlll99^ . We also ran models with 



the MESA code (jPaxton et al.ll201lL Samuel Jones, pri- 
vate communication) to study the lower mass limit for 
core-collapse supernovae and obtained similar results 
concerning the impact of the triple-a reaction rate on 
the boundary mass for carbon burning. Based on the 
models we have calculated, we estimate that the OKK 
rate decreases the minimum mass for core-collapse su- 
pernova by up to ~ 2Mq. 

Finally, we comment on the comparisons of our models 
with those computed bv lMorel et al.l (|2010f ). Our models 
of 4 - 6 Mp, stars are consistent with those of [Morel et ahl 
(J2010D . although the metallicity range is different. We 
also computed 6Mq models with solar composition to 
confirm the blue loops on the RGB. We find two blue 
loops during the central helium burning and during the 
growth of C -I- O core, though the size of the loop is 
smaller than their models. For [Fe/H] = —3, both blue 
loops are not observed because the helium burning starts 
during the contraction of the helium core on the subgiant 
branch, which prevent from the appearance of the first 
loop. The second loop appears when the mass of the C -I- 
O core is smaller than Schonberg-Chandrasekhar mass 
limit, i.e., w 0.37 times the mass of the helium core. For 
low metallicity stars, the C -I- O core mass exceeds this 
limit on the RGB, which results in the absence of another 
loop. If we adopt the OKK rate, our QMq model agrees 
qualitatively well with their 6.6Mq model with the OKK 
rate. 

3. EFFECT OF THE TRIPLE-a REACTION RATES ON 
HELIUM IGNITION 

In the previous section, we see significant modifica- 
tions to the current understanding of stellar evolution by 
adopting enormously enhanced triple-a reaction rates. 
Here we want to establish constraints on the triple-a re- 
action rate based on the occurrence of the FGB, which 
is supported by observations. We derive the maximum 
cross section for the triple-a reaction by considering its 
effects on the temperature of helium ignition, which cor- 
responds to the maximum luminosity of the FGB. In the 
following, we derive the helium ignition point from stel- 
lar evo l ution theory and model s based onlFuiimoto et a,l.l 
([T98l : [Fmmotol(ll982al[R l and lSuda fc Fuiimotol(l2010D . 
respectively, and estimate the allowed range of the triple- 
a reaction rates that is compatible with the observations. 

The ignition of hel ium shell burning is discussed by 
iFuiimoto et al.l (J1981D and its occurrence is described in 
eq. (9). Here we may simply define the helium ignition 
point by the requirement that helium ignition occurs in 
a shell where the nuclear timescale reduces to be as small 
as the thermal diffusion timescale. According to eq.(7) 
in IFuiimoto et al] (|1981D . the diffusion timescale is in in- 
verse proportion to the rate of diffusion, e^, which is given 

by 



Cr^ 



Aac T^ 

AacG M^T^ 
3 Kr4p2 



(1) 
(2) 



where Hp is the pressure scale height and the other sym- 
bols have their usual meanings. In deriving the condition 
for helium ignition, we assume that the rate of diffusion 
is identical to the rate of nuclear energy generation, e„. 
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at the shell potentially producing the largest nuclear en- 
ergy, i.e., the shell at the highest temperature. Therefore, 
we can obtain the ignition cross section given by the re- 
placement of e„ with er in the expression for the cross 
section of the triplc-a reaction. 



Ni{av)'^r 



1 



6 



NaYI p2 E^^^ ■ 



(3) 



where Y4 is the number abundance of ^He, and Eaaa the 
Q value of triple-a reacti on, which is set at 7 .274 MeV 
in this study according to lAngulo et al.l ()1999l ). 

The top panel in Figure |3] shows an example of the 
ignition curve for IMq model with Z = 0.02 using the 
standard set of NACRE reaction rates. Each point de- 
notes the helium ignition point as a function of the max- 
imum temperature in the model from the onset of hy- 
drogen shell burning to the beginning of the core helium 
flash. The maximum temperature occurs at the center 
until the core mass grows larger than Mi = 0.404Mo 
for logTmax = 7.80, as shown by a break in the lines. 
After that, it is shifted to the outer shell as the cen- 
tral region is cooled effectively via neutrino losses and 
develops a temperature inversion. The corresponding lu- 
minosity and helium core mass are shown in middle and 
bottom panels respectively. This model experiences a he- 
lium shell flash at the mass co-ordinate Mr = 0.230Afo 
when Ml = 0.477Mq. The density and temperature in 
the igniting shell are 10^-^^ g cm~^ and 10^'^'' K respec- 
tively. As shown in the figure, the ignition point is well 
predicted by the intersection of the ignition curve and 
nuclear reaction cross section provided by the NACRE 
compilation (dotted curve) as shown by the right vertical 
line in the figure. The prediction for the ignition point 
also holds for the OKK rate (dot-dashed curve). The 
derived ignition point for this model suggests that the 
helium ignition occurs at the temperature of the shell 
with T « 10^-^9 K if we adopt the OKK rate (see dotted 
vertical line in the left). The corresponding luminosity 
and core mass are log(L/L0) ~ 1.1 and Mi « O.2A/0 
respectively, shown in the middle and bottom panel of 
Figure |4] and emphasized by the vertical line. These val- 
ues are in quite good agr eement with the results of the 
same model computed bv lDotter &: Paxtoiil (|2009[ ). 

Figure [5] compares the ignition lines derived for the 
models of low-mass stars with various metallicity and 
metal-free stars. The different trajectories at different 
metallicities in the range 7.4 < logr(K) < 7.7 reflects 
the different thermal evolution of the core during the 
subgiant branch where the thermal energy in the hydro- 
gen burning shell is larger than the gravitational energy. 
This leads to the variations in the intersection of the cross 
sections between model predictions and the OKK rate 
since the mass of helium core differs during the transition 
from the central hydrogen-burning to the shell hydrogen- 
burning phase. Once the star evolves to the RGB, the 
rate of diffusion in the core is almost independent of to- 
tal mass and metallicity. Therefore, the values of the 
ignition cross section at the center are almost identical 
to each other at logT > 7.7 until the temperature inver- 
sion develops and the largest nuclear burning rate takes 
place in the off-center shells. A slight difference of the 
ignition cross section arises after neutrino energy losses 
become active, which simply results from the difference 



in the compression rate of the helium core due to hydro- 
gen shell burning. 

For 0.8Mq stars, the ignition temperature with the 
OKK rate becomes log T(K) w 7.5 - 7.6 with weak metal- 
licity dependence. The corresponding luminosity on 
the FGB in the H-R diagram can be seen in Figure El 
The comparison is made using the SAGA database 
(jSuda et al.ll2008() that contains the stellar parameters of 
low-mass halo stars with metallicity mainly in the range 
—4 < [Fe/H] < —1. The points on each evolutionary 
track show the helium ignition point at the temperature 
of logT(K) = 7.5 to 7.9 in 0.1 dex steps. In order to be 
compatible with the observations, the ignition tempera- 
ture has to be as large as logT(K) ~ 7.9. Accordingly, 
we find that the upper limit for the triple-a reaction cross 
section will be at most 2 to 3 orders of magnitude larger 
than the current rates as given by the NACRE compila- 
tion in order to remain coherent with obvsevations, un- 
less there are any other changes to the input physics. In 
other words, the shape of the cross section can be arbi- 
trary as a function of temperature as long as it intersects 
the ignition curve at the sufficiently high temperature to 
delay the helium ignition. More specifically, the cross sec- 
tion of the triple-a reaction should be less than ^ 10"^^ 
cm^ s~^ mole~2 at logT « 7.8. Another condition is 
that the cross section must have large temperature de- 
pendence of 1/ > 10 {{av) oc T'') at 8.0 < logT < 8.1 
at least to trigger the shell helium flashes as discussed in 
the previous section. These constraints arc emphasized 
by hatched areas in Fig. \E\ 

4. CAN OTHER UNCERTAINTIES COMPENSATE FOR THE 
CHANGE INTRODUCED BY THE OKK RATE? 

As discussed in the previous sections, astrophysical 
constraints clearly rule out the OKK rate. In this section, 
we try to work out if the effects of the very large OKK 
reaction rate can be balanced by changing other input 
physics and reaction rates. In particular, we focus on the 
reaction rates of the CNO cycle to preserve the FGB with 
variations in the triple-a reaction rate, such as the OKK 
rate. Other input physics such as opacities, equation of 
state, and neutrino loss rates are not worth considering 
because they will not contribute to dramatic changes of 
the evolution during the FGB. At least, there is no sup- 
porting evidence that stellar physics should be modified 
by more than an order of magnitude in the viewpoint 
of neither theory nor observations. For nuclear reaction 
rates, there are non-negligible uncertainties in some reac- 
tions, although the reactions related to hydrogen burning 
are rather well determined and the current rates used in 
stellar models arc compatible with observations. 

A possible solution to retain the FGB is the idea that 
enhancing the CNO cycle reaction rates accelerates evo- 
lution during the ascent of the RGB before igniting he- 
lium in the core. Since red giants are supported by hydro- 
gen shell burning mainly by CNO cycles, their enhanced 
reaction rates be able to increase the luminosity of the 
star rapidly. In order to estimate how much enhance- 
ment is required to accelerate the FGB evolution up to 
log L/Lq ~ 3, we compute models with enhanced reac- 
tion rates for the CNO cycles. The maximum luminosity 
on the FGB obtained by the model computations is plot- 
ted by filled circles in Fig [71 One should note that all 
the reaction rates related to CNO cycles have to be mul- 
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tiplicd by the same factor in order to increase the max- 
imum luminosity. Otherwise the bottle neck reactions 
prohibit the rapid growth of helium core and terminate 
the ascent on the FGB. In our models, the maximum 
luminosity saturates for the enhancement factors of 10^" 
or larger and cannot reach > lO'^i©. Therefore, there is 
no solution to attain log L/Lq = 3 as long as the OKK 
rate is used, even if we adopt enhanced CNO cycle re- 
action rates much larger than realistic values. For other 
metallicities, we estimate the maximum FGB luminosity 
using a semi-analytic method shown by lines in Fig. [71 
which is discussed in the Appendix. 

Figure H] shows the computations of 0.8Mq models 
with [Fe/H] = —3 by adopting the OKK rate with en- 
hanced CNO cycle reaction rates. In addition, our mod- 
els of low-mass stars with enhanced CNO reaction rates 
display quite different behaviors on the H-R diagram; 
the turn-off point at the main sequence moves toward 
the right because of the dominance of CNO cycle at the 
beginning of the main sequence. 

5. SUMMARY AND CONCLUSION 

We investigated the impact of changing the reaction 
rate of the triple-a reaction on stellar structure and 
evolution. Models with the reaction rate provided by 
I Ogata et all ([200^^ ) (the OKK rate) are computed for 
stars in the mass range from 0.8 to 25Mo. It is shown 
that the OKK rate brings about several serious problems 
in astrophysical contexts including the disappearance of 
th e first red giant branch (FGB) as already pointed out 
bv lDotter fc Paxtoiil (J2009D . Here we have found another 
significant problem that low-mass stars do not experi- 
ence helium shell flashes due to the weak temperature 
dependence of the reaction rate. The absence of helium 
shell flashes results in the absence of the change of sur- 
face chemical composition through the third dredge-ups. 
This would suppress important channels to the forma- 
tion of carbon-enhanced stars since, for a majority of 
carbon stars, their abundances are best explained by 
AGB evolution or by binary mass transfer from AGB 
stars. In addition, it would also prevent the formation 
and transportation mechanisms of ,s-process elements in 
low- and intermediate-mass stars. As a consequence, 
AGB stars would no longer be the progenitors of carbon 
stars and ,s-process element enhanced stars. In partic- 
ular for metal-poor halo stars, the large fraction of ob- 
served carbon-enhanced stars with s-process element en- 
hancement could not be explained by mass transfer from 
the former AGB stars. Observationally, the detection of 
Technetium in the surface of AGB stars is the support- 
ing evidence of the dredge-up of s-process elements by 
the deepening of the convective envelope. On the other 
hand, for massive stars with M > IOMq, the effect on 
the evolution is small because the helium ignition tem- 



perature is so high tha t the difference in th e reaction rate 
between the NACRE (jAngulo et al.lll999f ) and the OKK 
rate is not significant. 

We also quantitatively estimated the allowed range of 
the triple-a cross section under the currently adopted 
input physics in stellar evolution code. In our analysis 
of the helium ignition point of the star, it is suggested 
that the triple-a reaction rate should not be larger than 
~ 10~^^ cm^ s~^ mole~^ at 10^'^ K in order for stars 
to still have a red giant phase. In addition, it is shown 
that the occurrence of helium shell flashes requires the 
temperature dependence of the triple-a reaction cross 
section with at least z/>10atT~l-1.2xlO*K where 
the cross section is proportional to T^. For other tem- 
perature ranges, the shape of the reaction rate can be ar- 
bitrary as a function of temperature only if it intersects 
the ignition cross section given in eq.([3|) at logT > 7.9 
where helium ignition occurs at the FGB luminosity of 
L> IO^Lq. 

A way to compensate for the changes due to the OKK 
rate was sought by pursuing the increase of the CNO cy- 
cle reaction rates, but found to be impossible. We looked 
for a solution to increase the FGB tip luminosity by in- 
creasing the CNO cycle reaction rates to accelerate the 
growth of helium core. Even if we adopt the enhance- 
ment factor of > 10^" for all the reaction rates related to 
CNO cycles, the maximum luminosity of the FGB cannot 
reach lO^i©. 

The met hod to derive t he trip le-a reaction rate pro- 
pos ed bv iNornoto et al.l (J1985I ) is severely criticized 
bv lOgata et al.l (I2009D . while the result obtained by 
lOgata et al.l ( 2009() is still far from being consistent with 
observations. Since the triple-a reaction is in the wide 
range of applications such as accretion events on com- 
pact stars like novae. X-ray bursts and the determination 
of cosmological distances through the Cepheid variables, 
more investigations for this important reaction arc de- 
sired. 
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APPENDIX 



SEMI-ANALYTIC METHOD TO DERIVE THE MAXIMUM LUMINOSITY OF THE FGB 

As a semi-analytic test to quantify the maximum luminosity at the tip of the FGB, we derive the relation between 
CNO cycle reaction rates and the maximum luminosity. Using the expres sion for nuclear energy generation rate by 
e„ oc p'T"^, the luminosity is given by the following thin-shell approximation (jHavashi et al.ll963 : lSugimoto fc Fuiimotol 
|1978() using the homologous invariant, U and V, defined by f/ = 'iirr^p/Mr and V = GMrp/{rP), respectively 
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where subscript 1 denotes the hydrogen-burning shell and ejv the nuclear energy generation rate by CNO cycles, 
total stellar luminosity is approximated by the quantity at the hydrogen burning shell, which is written as 
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where /3i is the ratio of gas pressure to total pressure at the edge of the helium core. By equating eq. (IA2[) and (jASp . 
together with the assumption that Vi = 4 and the equation of state, we obtain a set for solutions of density, temperature, 
and pressure for a given mass and radius of the helium core, and che mical composition. The assumption that Vi = 4 a t 
the hydrogen burning shell is based on the theory of shell burning (jHavashi et al.|[r962l : iSugimoto fc Fuiimotol [19781) . 
We also assumed the polytropic index of A'^ = 3 and the surface chemical composition as the input abundances in 
solving the equation of state. The mass. Mi, and radius, ri of helium core at the ignition of helium burning are taken 
from the model derived for the OKK reaction rates. Then, the maximum luminosity, L^q^^-, reached in the FGB 
with the OKK reaction rate is obtained as a function of enhanced factor, /cnOi of CNO cycle nuclear reaction rates 
by the above procedure, which is shown in Figure [T] For a given mass and radius of the helium core, the increase in 
the luminosity stems solely from the contribution, 1 — /3i, of the radiation pressure to the total pressure, as seen from 
eq. (jASp . The increase in the dominance of the radiation pressure is realized through the decrease of the density in the 
burning shell because of the strong temperature dependence of the nuclear burning rate. Accordingly, the dependence 
of the maximum luminosity on the enhancement factor of CNO cycle is weak, as shown in Figure [T] 

The prediction obtained by the semi-analytic model shows a good agreement with model computations for the 
enhancement factor up to /cno — 10^ with the logLpQ^^- < 2.5. For still larger enhancement factor, however, the 
numerical models no longer give an increase in the maximum luminosity on the FGB, differently from the analytical 
predictions. The decrease of maximum luminosity in numerical models arises from the expansion of the helium core 
caused by the decrease in the density in the hydrogen burning shell with a finite temperature, while the value of ri is 
fixed in the above analytical estimate. 
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Fig. 1. — Comparison of selected evolutionary tracks on the H-R diagram using the OKK (solid lines) and NACRE (dashed lines) rate. 
Initial mass (in Mq) and metallicity (in [Fc/H]) of stars arc labeled on the next to the main sequence. For 12 and 25AIq with [Fc/H] = —0.3, 
we plot the models computed with the Geneva code. 
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Fig. 2. — Comparison of helium burning luminosity as a function of helium core mass for the models with the NACRE rate (top panel) and 
the OKK rate (bottom panel). The definition of the lines is the same as in Fig. [T] but the models undergoing carbon burning (M > 6Mq 
with the OKK rate) arc excluded. 
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Fig. 3. — Cross sections for triplc-o reactions and their temperature dependences {d\nl^a v) aaa I d\\iT , see right axis) for the NACRE 
rate and the OKK rate. 
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Fig. 4. — Model result of IMq with Z = 0.02. Top panel shows the compariso n of the exp e cted ignition points computed by eq. Jsjl 
with the cross sections by triple-a reaction taken from NACRE (broken line) and lOgata et all II2009I . dash-dotted line), as a function of 
temperature at the shell where the triple-a reaction reaction rate is maximum. Middle and bottom panel are the luminosity and helium 
core mass, respectively. Dotted vertical lines represent the intersection between the reaction cross sections and ignition curves determined 
from stellar models. 
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Fig. 5. — Ignition curves derived by our stellar models compared with the cross section of triplo-o reaction of NACRE and lOgata et aO 
l|2009l ). Hatched areas denote the constraints on the reaction rate obtained by this study. Hatched circle designates the maximum cross 
section to avoid central helium ignition in lowr-mass red giants. Hatched square showrs the temperature range for which the temperature 
dependence of the cross section should be large enough {u > 10) so that low- and intermediate-mass stars undergo helium shell flashes. 
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Fig. 6. — The end point of the FGB for models with the different triplc-o reaction rate on the H-R diagram. Each filled circle on 
evolutionary tracks shows the ignition point when the ignition temperatur e is logT = 7.5, 7 .6, 7.7, 7.8, and 7.9 from bottom to top on the 
same line. Observed halo stars are also plotted using the SAGA database l|Suda et al.ll2008l ). 
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Fig. 7. — The maximum luminosity in the FGB obtained by multiplying the cross sections of CNO cycles by the factor /cNG- Each 
circle denotes the luminosity at the tip of the FGB obtained by the computation of a model of O.8J\/0 with [Fe/H] = —3 using the OKK 
rate with enhanced CNO cycle reaction rates. Linos show the predictions by our semi-analytic model (see the Appendix). 
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Fig. 8. — Evolution of O.SMq stars and [Fc/H] = —3 with the OKK rate and enhanced CNO cycle reaction rates. The fiducial model 
and the models with OKK rate and without the enhancement of CNO cycle reaction rates arc plotted for comparison{/cNO = !)■ The 
evolution is followed until the helium core grows nearly to the total stellar mass for the models with enhanced CNO cycle reaction rates. 



